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Abstract 

The recently introduced deflected mirage mediation (DMM) model is a string-motivated par- 
adigm in which all three of the major supersymmetry-breaking transmission mechanisms are 
operative. We begin a systematic exploration of the parameter space of this rich model context, 
paying special attention to the pattern of gaugino masses which arise. In this work we focus 
on the dark matter phenomenology of the DMM model as such signals are the least influenced 
by the model-dependent scalar masses. We flnd that a large portion of the parameter space in 
which the three mediation mechanisms have a similar effective mass scale of 1 TeV or less will 
be probed by future direct and indirect detection experiments. Distinguishing deflected mirage 
mediation from the mirage model without gauge mediation will prove difficult without collider 
input, though we indicate how gamma ray signals may provide an opportunity for distinguishing 
between the two paradigms. 

1 Introduction 

If supersymmetry is relevant at the electroweak scale then forthcoming experimental data should 
reveal the presence of new states which will be studied extensively at colliders. If the lightest 
supersymmetric particle (LSP) is stable it will provide an excellent candidate for explaining the 
presence of non-baryonic dark matter in the cosmos [1] . The phenomenology of these new states - 
both at colliders and in cosmological observations - will be determined by their interactions with 
each other and with the states of the Standard Model. These are determined, in turn, by the 
manner in which supersymmetry breaking is communicated from a hidden sector to the observable 
sector [2]. 

Over the past ten years most studies of supersymmetric phenomenology have tended to focus on 
one of three general paradigms for mediation of supersymmetry breaking to the observable sector. 
These are modulus (or "gravity" ) mediation [3l S] , gauge mediation [5] and anomaly mediation [6l [7] . 
The latter two are often motivated by bottom-up concerns such as flavor-changing neutral currents 
and CP-violation, but the case of modulus mediation (a special string-motivated case of general 
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supergravity mediation) can avoid phenomenological constraints as well, provided the modulus or 
moduli in question have family-universal couplings to the observable sector and the underlying 
string theory has certain well-motivated isometries [SI [S]. These paradigms are typically studied 
in isolation, assuming that only one such mechanism dominates the pattern of soft supersymmetry 
breaking. As an example one could consider the hugely influential "Snowmass Points and Slopes" 
benchmark models designed for studies of collider phenomenology fTO| . 

This is a reasonable starting point, but recent work suggests that there are well-motivated rea- 
sons to consider pairs of these mediation mechanisms - or perhaps all three mechanisms - simulta- 
neously. The combination of anomaly mediation and gauge mediation, known as deflected anomaly 
mediation, was studied soon after the anomaly-mediated scenario was first proposed [111 [T2] . For 
the combination of anomaly mediation and modulus mediation to produce comparable contribu- 
tions to soft supersymmetry breaking a hierarchy of mass scales must be engineered. This was first 
observed in certain classes of heterotic orbifold models |13| [TH [15] but received increased attention 
with the advent of KKLT-type moduli stabilization in D-brane models |16| [T7 j. A particular class 
of Type IIB string compactifications with fluxes [TH] gave rise to the so-called "mirage mediation" 
models [H |20l [M] whose phenomenology has been extensively studied [22l l23l IM} [25] [26] [271 128] . 

Recently the two combinations described above were combined in a scheme dubbed "deflected 
mirage mediation" (DMM) in which all three paradigmatic transmission mechanisms are operative 
and comparable [29] While the DMM framework does not result in fully generalized superpartner 
mass possibilities - particulary for the all-important gaugino sector, as we will see below ~ it does 
provide a very rich phenomenology. Like the mirage mediation case which preceded it, some basic 
elements of the phenomenology of DMM have been studied within the contexts of specific string- 
based frameworks |30| \3T\ [32] . In this paper we embark on a program to study the following 
general question: if all three supersymmetry-breaking mediation mechanism are present in roughly 
equal amounts, how will experimental observation reveal this fact? By phrasing the question in 
this manner we hope to study how certain "underlying principles" (such as the mechanism of 
supersymmetry-breaking mediation) can be extracted directly from data without reconstructing 
the full low-energy Lagrangian |34| [35] . 

Here we therefore wish to study primarily the gaugino sector, which reflects fewer model- 
dependent properties. Indeed, the DMM scenario continues to be a two-parameter family of gaugino 
mass patterns, with one parameter representing deviations from universality and the other the over- 
all mass scale. In this paper we will consider the dark matter signals of the deflected mirage media- 

^See also the construction in '33'. 
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tion scenario, which are determined largely by the physics of the gaugino sector alone. Of course full 
model-independence is impossible. Such important properties as the masses of the supersymmet- 
ric Higgs fields and the Higgsino content of the LSP will depend on soft supersymmetry-breaking 
scalar masses which are always model-dependent. We will consider the truly model-dependent 
implications for hadron collider physics in a future publication [36] . 

In its simplest manifestation the general mediation scenario would involve three mass scales, one 
each for the (bulk) gravity contribution, superconformal anomaly contribution and the gauge mes- 
senger contribution. This can be re-parameterized by a single overall scale and two dimensionless 
ratios. Such was the parameter set adopted by Everett et al. [29]. We will review the construction 
of the deflected mirage mediation model in Section [2] focusing most strongly on the gaugino sector 
where the interplay of supersymmetry-breaking mediation mechanisms is most transparent. We 
will study how the addition of gauge mediation alters the allowed parameter space of the theory. 
We then turn to dark matter observations in Section [3j indicating how the deflection arising from 
gauge-charged messengers alters the phenomenology relative to the un-deflected mirage model. Our 
survey will include a number of current and future experiments seeking to detect the presence of 
relic neutralinos. We will find that regions of the parameter space which result in a wino-like or 
mixed Higgsino/gaugino LSP will be probed in direct detection experiments at the one ton- year 
scale as well as in a variety of indirect detection experiments, provided that the effective mass scales 
of the various mediation mechanisms are in the 1 TeV range or less. These are also the areas in 
which the thermal relic abundance of the LSP is no greater than the upper bound inferred from the 
WMAP data. Regions of the parameter space with a heavier spectrum, or with a predominantly 
bino-like LSP, will prove much more difficult to observe via cosmological or astroparticle observa- 
tions. We perform a basic survey of cosmic ray signals for the DMM model and find that none of 
the parameter space gives rise to an adequate explanation for the recent ATIC data, though some 
points can give a reasonable fit to the PAMELA data with boost factors less than 0(100). How- 
ever, applying such boosts to the anti-proton flux as well would make these points in conflict with 
a lack of signal in the anti-proton data without modification of the diffusion model. In Section |4] 
we consider the extent to which the observations discussed in Section [3] will be helpful in distin- 
guishing the deflected mirage mediation model from the previously-studied mirage pattern without 
gauge-charged messengers. Here we will encounter a "dark matter inverse problem" in trying to 
distinguish between the two cases. We will demonstrate why this occurs and suggest some possible 
avenues for future study to resolve this problem. 
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2 Parameterizing Deflected Mirage Mediation 



In this section we will review the construction of the soft Lagrangian in deflected mirage mediation. 
As this material has been presented more fully elsewhere [3^ [32] we will be brief in our discussion 



of the underlying string-theoretical context. We begin in Section 2.1 with a basic derivation of the 
gaugino mass patterns in deflected mirage mediation in a model-independent manner. This will 
allow us to introduce some important notation and outline the parameterizations we will be using 



throughout the rest of the work in Section 2.2 In Section 2.3 we complete the construction by 
introducing the soft supersymmetry-breaking scalar masses and trilinear A-terms and explore the 
allowed regions in the parameter space. 



2.1 Gaugino Masses in Terms of Mass Scales 

In this subsection we present the construction of the gaugino mass patterns without any reference 
to its possible origin from string-theoretic considerations. The treatment here closely parallels that 
of [35] in the absence of gauge-charged messengers. For the sake of theoretical clarity we will work 
to one-loop order in the renormalization group equations in this section. We begin by assuming 
three contributions to the soft supersymmetry breaking gaugino masses of the MSSM. Let the 
contribution from Planck-suppressed operators be universal in form with a mass scale given by the 
quantity Mq. The contributions from the superconformal anomaly will be proportional to some 
(a priori different) mass scale given by Mg. We assume that these contributions arise at some 
effective high-energy scale /Xuv at which supersymmetry breaking is transmitted from some hidden 
sector to the observable sector. It is common to take this scale to be the GUT scale at which gauge 
couplings unify, but in string constructions one might choose a different (possibly higher scale) at 
which the supergravity approximation for the effective Lagrangian becomes valid. Finally, we will 
include a contribution from some gauge-charged messenger sector which is proportional to a third 
mass scale Amess- As in the original work of |29| |30] we will assume this messenger sector comes in 
complete GUT representations of the Standard Model gauge group so as to preserve the successful 
gauge coupling unification of the MSSM. In particular we assume Nm copies of 5,5 representations 
under SU{5), which give rise to contributions to soft supersymmetry breaking gaugino masses at 
the energy scale /^mess ^ Muv- 

The full gaugino mass at the high energy boundary condition scale /iuv is given by the expres- 
sion [la ET] 

b' 

Ma (/xuv) = Mo + gl M jf^m3/2 , (1) 
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where a = 1,2,3 labels the Standard Model gauge group factors Qa and is the value of the 
gravitino mass, given by = (e^^'^W\. The quantities b'^ are the beta-function coefficients for 



the Standard Model gauge groups. In our conventions these are given by 

ba = -{3Ca-J2C:), (2) 

i 

where Ca, C* are the quadratic Casimir operators for the gauge group Ga, respectively, in the 
adjoint representation and in the representation of the matter fields charged under that group. 
As the summation in Q is over all degrees of freedom in the theory present at the scale /iuv > /^mcss) 
it must include the contribution from the gauge-charged messenger sector. Therefore we may write 

b'^ = ba + Nm (3) 

with ba the beta-function coefficients of the MSSM in the absence of these messenger fields 

{&l,&2,63} = |f ,1,-3| . (4) 

Note that in the absence of the messenger sector {Nm = 0) if we take /iyv = /^gut then we would 
have 

(/^uv) = 92 (muv) = gi ifj-vv) = glvT (5) 

In the presence of the gauge-charged messengers the unification of gauge couplings continues to 
occur at a scale //gut — 2 x 10^^ GeV but the value of the unified gauge coupling differs from that 
of ([5]) by the relation 

1 1 A^™^ //^\ 



'rness 



5a (^uv) 5gut 87r2 \n 
which is valid at one-loop order for all values of a. 

We now imagine evolving the expressions in ([T]) to the scale //mess via the one-loop renormal- 
ization group equations. The first term in ([T]) is renormalized by a multiplicative factor at one 
loop 



M*--^ (//mess) = Mo 



b ( 

1 - gl (^mcss) ^ In 



//mess 



(7) 



while the second term in ([T]) can be evolved by simply replacing the gauge coupling with its value 
at the intermediate mass scale 

M^^^^ (^^^^^^ ^ ^2 J^m3/2 . (8) 
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Combining ([7| and Q we have 

Ma (Mmess) = ol (/^n 



6' 



167r^ 



3/2 



+ Mo 



1 - 5a (/^mcss) 



87r2 



In 



/^mess 



(9) 



Clearly, if the two parts of the expression in ([9| are to be of roughly the same size it will be 
necessary to engineer a situation in which m-^i2 — IQtt'^Mq. 

At this intermediate scale we may now integrate out the gauge-charged messenger sector, which 
in the presence of supersymmetry breaking parameterized by m3/2 generates a threshold correction 
to the gaugino masses of the form |38j 



AMa 



167r^ 



+ m 



3/2 



Evolving the combination Ma (/imess) + from the intermediate scale /in 

scale /Xew (again at one-loop) results in the final expression 



(10) 

to the low-energy 



Ma{ll^ 



al (/^: 



Mo ( 1 



m 



3/2 



+ 



1 - al ifJ-Ew) 



ha 



In 



/^mess 



alif^B 



167r2 



-A, 



From (11) we immediately see that we must engineer the relations Aj, 



"i3/2 



(11) 



IGvr^Mo if ah 



three contributions to the gaugino masses are to be comparable in magnitude. Let us note before 
going forward that the preceding derivation made tacit assumptions about the relative phases 



between the various contributions in the final expression (11). Given a true top-down construction 



these phases should be determined by the same mechanism which stabilizes the auxiliary fields in 
the supergravity Lagrangian. 

2.2 A Theory-Guided Parametrization 



The masses in (11) depend explicitly on the three mass scales Mq, 7713/2 and Amess as well as the 
number of messenger multiplets N^- They also depend implicity on the value of the intermediate 
scale /imess through the logarithms and gauge couplings. This somewhat unwieldy expression can 
be simplified by relating the explicit mass scales to some overall mass scale via two dimensionless 
ratios. Following Everett et al. we can define 



Ar 



m 



3/2 



(12) 
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which will naturally be of order unity for the cases of interest to us here. To obtain another ratio 
of 0(1) it is common to define 

Oim = 7 — V , (13) 



Moln(^MpL/m3/2j 

where MpL is the reduced Planck mass MpL = 2.4 x lO"*^^ GeV. Again, for the models of interest to 
us here we will have Um typically of order unity. Using the definitions (12) and (13) the expression 



in (11) can be re-expressed as 



Mo 



1 - al (a^ii 



87r2 



In 



A'mese 



1 - gl (Mew) ^ In 

(/^mess) 



• (14) 



With these conventions one can define the so-called "mirage scale" /Xmir at which the three soft 
supersymmetry breaking gaugino masses have an identical value as 



Mmir — /^GUT 



m 



■3/2 



Mp 



where the parameter p is given by 
P 



^ _^ Nmgl (A^Uv) / fJ'GVT 
8vr2 VA^mcss 



/ 



^ _ amagNmgl (/^ 



167r^ 



m 



3/2 / _ 



(15) 



(16) 



The expression in ( 14 ) is still somewhat cumbersome. It was noted by Choi [31j that the gaugino 



mass expressions (11) at the low scale /Uew can be re-arranged into the form 



Afa(AiEw) = 

where we have introduced the variable 



/^mir 



/3a(/i) 



and have defined the effective mass scale 



M, 



eff 



RMq ; R = l 



87r2 



87r2 



Mpi 



+ In 



m 



3/2, 



A'gut 
Mmess 



(17) 



(18) 



(19) 



The expression in ( 17 ) is precisely identical to that of ( 14 ) at one loop order, provided we remember 
that the high-scale gauge coupling g^^jrj, appearing in (19) is related to the high-scale gauge coupling 



7 



9a (Muv) appearing in (14) via ([6|). We note that the dimensionless quantity R is simply the inverse 



of the parameter p in (16), and thus 



/^GUT 



m 



3/2 



am/2R 



(20) 



In [32] the expression in (17) was further modified by introducing two new dimensionless vari- 

(21) 



ables based on the quantity i? in (19) 



1 



Oir. 



R + ar. 



R + Ur, 



which can be used to write (17) in the form 



Ma{^l) = Mo (1 + l3a{lJi)t) |l + y 
where we have defined two scaling variables t and t' via 



1 + PaW 



1 



(22) 



In 



/^GUT 



t' 




(23) 



The expression in (22) is intended to be evaluated at the energy scale /i = /Uew < Mmess- The free 



parameters in (22) are {Mq, x, y}, though (22) is still ultimately a function of only one scale and 
one dimensionless variable (as we will explicitly show in Section |4]). Having specified x and y values 
one may obtain Um = y/x, and then the value of m^i2 can be obtained from that of Mq by finding 



the solution to the definitional equation for am in (13). Though the series of definitions in (19) 



and (21) may appear to sacrifice clarity for the sake of brevity, they allow an extremely powerful 



way to survey all possible mediation mechanisms in a simple two-dimensional plane. 



2.3 Exploring the Parameter Space 

As we have seen, the gaugino mass sector in deflected mirage mediation is particularly transpar- 
ent. Though the three soft supersymmetry- breaking gaugino masses are sufficient to determine a 
large part of the dark matter phenomenology of the LSP neutralino, the scalar sector will also con- 
tribute in important, indirect ways. These include determining the /i-parameter via the electroweak 
symmetry-breaking (EWSB) conditions and thereby the LSP wave- function, and by allowing cer- 
tain important resonant annihilation or co- annihilation processes in the early universe. Therefore 
we cannot adequately discuss the prospects for dark matter signals in deflected mirage mediation 
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without specifying the rest of the soft supersymmetry breaking Lagrangian, particularly the scalar 
masses and trilinear couplings. 

Doing so requires working within a model framework. We will here very briefly outline that 
framework and the relevant expressions needed to go forward. Further background can be found 
in the relevant references, particularly |30| l32] . We imagine Type IIB string theory compactified 
on a Calabi-Yau orientifold. To this is added flux in the manner of KKLT, producing a region of 
highly warped geometry. The remainder of the compact space remains (approximately) a Calabi- 
Yau manifold. These background fluxes serve two purposes: the first is to stabilize most of the 
geometrical moduli associated with the construction. The second (and related) purpose is to 
produce a constant in the superpotential for the remaining modulus of sufficient size to produce 
stabilization with ~ 0(1)- This remaining modulus is a Kahler modulus which can be stabilized 
by field-theoretic mechanisms such as gaugino condensation. The observable sector {i.e. the MSSM) 
resides on Dj, or D-j branes (or both) in the bulk of the internal manifold. In this minimal set-up it 
is plausible that a conspiracy of scales can be achieved in which direct modulus-mediation (via non- 
vanishing F-terms for the Kahler modulus) and loop-induced anomaly-mediation are competitive 
contributions to observable sector scalar masses. 

Let us denote the (chiral) Kahler modulus superfield by T and its highest auxiliary component 
by Ft- If the gauge kinetic functions fa for the three Standard Model gauge groups are universally 
given by fa = T, then a non-vanishing vacuum value for Ft will generate a universal component 
for the gaugino masses as in ([T]) 



where t = T\q=q is the lowest (scalar) component of the superfield T. The quantity in ( |24[ ) will 
also generate scalar masses and trilinear scalar couplings. The exact form depends on how the 
Kahler modulus appears in the Kahler metric for the matter superfields, if at all. Such couplings 
are difficult to compute, but their leading behavior at large volume (large {t + -£)) can be readily 
inferred and parameterized. We will assume a Kahler potential of the form 



with modular weight rii for the matter superfield ^i. These numbers depend on the location of the 
observable sector matter fields and may be vanishing. It was argued in j39| , using the work of , 
that for intersecting D-brane models in the Type IIB context one should expect nj = if $j is 
localized on branes, = 1 if $j is localized on D^ branes, and rij = 1/2 for matter localized at 




(24) 



K = ^(r-hr)-"^$i$ 



1 1 



(25) 
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the intersection of branes. We point out, however, that no exphcit model realizing both the fluxed 
stabilization mechanism and the MSSM matter content has yet been constructed. 



In the absence of gauge-charged messengers, the contributions of (24 ) and the assumption in (25 ) 



are sufficient to specify the form of the scalar sector soft supersymmetry breaking terms at a scale 
Huy. Such a system defines the mirage mediation (or mixed modulus /anomaly mediation) scenario. 
To this we now wish to add Nm copies of gauge-charged messengers in the representation 5,5 
under SU{5). The presence of such terms introduces new contributions to the soft scalar masses 
and gaugino masses at the scale Hmess at which these fields are integrated out of the spectrum. 
The calculation of such terms follows the general procedure of gauge mediation models [5]. Soft 
supersymmetry breaking is transmitted to the observable sector states via loop diagrams involving 
the messengers. The size of these contributions is determined by the parameter Amess appearing 



in ( 10 ). In principle this mass scale can be of any size as it is a priori independent of the physics of 
moduli stabilization. The relative size of Amcss to (and hence the parameter ag) will depend 
on the model one postulates for dynamically generating masses and splittings for the messenger 
sector. Circumstances for which ag is of order unity are not hard to construct |29| [5^1 . 

The soft terms for the scalar masses and trilinear couplings, in the presence of all three 
supersymmetry-breaking transmission mechanisms, are given set the initia.1 scale /Xuv ^ Mmess 

by [30] 



(3 - n, - 
(1 - rii) 



nk) 



{li + 7i + Ik) 
167r2 



Omln (^MpL/m3/2 



In (^MpL/mg/aj 



7. 

(16^2)2 



In (^MpL/mg/a)) 



(26) 
,(27) 



with ara being the parameter introduced in (13) and the quantities 7i, 0[ and 7^ being functions of 



Standard Model gauge and Yukawa couplings. Their definitions and explicit values can be found 



in ^30^. The primes on the quantities in (27) indicate that these parameters are computed taking 
into account the presence of the messenger sector as well as the MSSM matter content. At the scale 
A* = Aimcss the messenger sector is integrated out of the theory, producing the contribution (10) for 
the gaugino masses 



AMa = -N^Mq- 



167r2 



(1 + ag) In (MpL/mg/a) , 



and a contribution 



Am? = M2^(2iV^C7,) 



(167r' 



2^2 



am(l + ag) In (^Mpi^/ 7113/2 



(28) 



(29) 
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for the scalar masses. These expressions can be re-expressed in terms of the parameter R of (19) 



and the set {x, y} of (21), but doing so is not particularly illuminating. 

Unlike the case of gaugino masses, no analytic expression at low scales is available for the scalar 
masses and trilinear couplings - even at one loop order - except in certain approximations [32]. 
Therefore when we wish to consider the full model-dependent set of soft parameters we must connect 
energy scales using the renormalization group (RG) equations. We do so in two stages. The first 
connects /Liuv = /^gut and /imess with boundary conditions given by Q, (|26j) and (|27j). At this scale 



the corrections (28 1 and (29) are added and the soft terms are evolved from Hmess to n^w = Mz- 
From here we pass the low-scale soft parameters to SuSpect 2 . 4 [l2] to check for proper electroweak 
symmetry breaking and to compute the physical masses of the various superpartners. 

We will investigate the parameter space defined by the set {x, y, Mq} by considering slices 
of the {x, y} plane for various values of Mq. We do this in two different regimes. In the first 
case we investigate the explicit model described in this section for the case where tan/3 = 10, 
Aimess = 10^*^ GeV, Nm = 3 and where the modular weights of the various MSSM states are given 
by 

{UQ, nu, UD, UL, UE, UH^, uhJ = {1/2, 1/2, 1/2, 1/2, 1/2, 1, 1} . (30) 

This choice of modular weights is not motivated by any particular construction; we choose it to 
match with certain benchmark models from the literature. We will call this the "model-dependent 
scenario." In the second case we will dispense with the model-dependent soft-terms by simply 
setting all trilinear couplings to zero, set ^ = rriA = 1 TeV and set all scalar masses to the value 
of the gluino soft mass M3 or to 1 TeV, whichever is larger. We do this by hand at the low energy 
scale. While artificial, this provides some sense of how the phenomenology of the gaugino sector 
alone is influenced by the choice of {x, y} for a given value of Mq. We will call this second case the 
"model-independent scenario." 

In Figure [1] we display the physically allowed region in the {x, y} plane for various values of the 



mass scale parameter Mq in the model-dependent case with modular weights given by (30). The 
primary constraint in the {x, y} plane comes from the direct search limits for gauginos. The mass 
bounds arising from these searches are somewhat model-dependent, particularly for the gluino. We 
will require that the lightest neutralino have a mass > 46 GeV and that the lightest chargino 
mass obeys > 103 GeV [l3]. For the gluino we will be much more conservative. Though 
recent searches at the Tevatron have become increasing model-independent, we do not wish to rule 
out too much of the parameter space of this model class prematurely. We will therefore require 
only that m- > 200 GeV [44J in our figures, though a bound of m~ > 300 GeV |35] will be more 
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500 GeV 




Figure 1: Allowed Parameter Space for the "Model-Dependent Scenario." Contours indicate the locus of 
points for which m~ — 46 GeV (dashed contour) , m~ = 103 GeV (solid contour) and m~ = 200 GeV (dash-dotted 
contour). The dark shaded region is the area in which the stau is the LSP, while the smaller hatched region in the 



center has a stop LSP for the modular weight choice of (30 1. The hatched region in the upper left of each plot 



indicates where no EWSB occurs. For larger values of the parameter Mo we have indicated the area in which the 
gluino (or the chargino) is the LSP by the darker shaded region(s). The labeled points are the benchmark models of 
Table [l] The intersecting lines indicate those points for which a™ = 1 (lower left to upper right) and where R — 
(upper left to lower right). The intersection of these two curves designates the prediction of the simplest KKLT 
scenario. 12 
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Model D 


Model E 


Model F 


Mo 


2000 GeV 


1000 GeV 


1000 GeV 


500 GeV 


500 GeV 


933 GeV 


"T.3/2 


62.6 TeV 


32.0 TeV 


32.0 TeV 


13.2 TeV 


51.3 TeV 


31.7 TeV 


/^mess 


10^2 GeV 


10^ GeV 


10^2 GeV 


lO^o GeV 


10^0 GeV 


N/A 


Om 


1 


1 


1 


0.81 


3.26 
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«9 


1 


-1/2 


-1 


0.14 


1.47 
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3 


3 


3 







0.668 


0.561 


0.472 


0.674 


0.401 


0.5 


y 


0.668 


0.561 


0.472 


0.543 


1.307 


0.5 



Table 1: Benchmark Models. The relevant mass scales are given in the first three entries, with the parameters 
of Everett et al. in the second block. The final block re-casts these parameters in terms of the parameterizations of 
Choi et al. from (21 1. All models are defined with positive value of /i and tan/3 — 10. Models A-B were considered 
in [29] while Model C was considered in [301. Model F is a mirage model without messengers near the prediction for 
the basic KKLT model. 



applicable over many regions of the parameters. The lightly shaded region is ruled out by these 
direct experimental constraints. In addition there is some parameter space for which the gluino 
(and for very large Mq, the chargino) can be the LSP, indicated by the darker shading. For this 
choice of modular weights the stau is the LSP for smaller values of x and y, indicated by the dark 
shaded area. There is also a small region where the lightest stop can be the LSP, indicated by the 
hatched region near the center of each figure. Finally in the upper left corner of each panel there 
is a region for which < for the pseudoscalar Higgs and no electroweak symmetry breaking 
occurs. 

This choice of parameter space offers a distinct advantage: the various theoretical model limits 
can easily be defined in terms of the dimensionless quantities x and y. In the {x, y} plane cases 
for which ii = lie along the line for which x + y = 1. The line extending from the origin defined 
by x = y is the case with am = 1 for arbitrary ag. The intersection of these two lines is the 
simplest prediction of the KKLT framework which inspired the mirage meditation scenario. We 
have imposed these theoretical model lines on top of the allowed space in Figure [T] The gauge- 
mediated hmit in Figure [l] is formally at the origin of the {x, y} plane where — > 0, but this 
limit must be taken while keeping the product R Mq finite |32j and is therefore impossible to reach 
in our figure. So too, the pure anomaly-mediated limit is impossible to reach since this is formally 
the point {0, 1} {i.e. Om — > oo) where the product Mq is held fixed. The pure modulus (or 
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Figure 2: LSP Properties for the "Model-Dependent Scenario." Contours of constant LSP mass are 



given for the scenario witli modular weights of (301 for mass scale A/q = 500 GeV (left panel) and Mo = 1000 GeV 
(right panel). The lightly shaded region is ruled out for the reasons indicated in Figure [l] For both values of Afo the 
LSP is primarily bino-like for y ^ 0.6 and y ^ 1.4. For the case of Mo — 500 GeV there is also some parameter space 
where the LSP is a mixture of Higgsino and wino, indicated by the hatched region. 



gravity) mediated limit is the point {1, 0} in the figures. For large enough values of Mq there 
is an expansive parameter space beyond the mirage mediation "frontier" for which a compressed 
spectrum of gauginos occurs. The marked points in Figure [T] represent specific parameter choices 
which we will single out for particular study later in the work. The parameter sets which define 
these cases are given in Table [TJ 

In Figure [2] we concentrate on the properties of the lightest neutralino for the specific cases of 
Mo = 500 GeV (left panel) and Mq = 1000 GeV (right panel). The union of all the theoretically 
and experimentally forbidden regions in Figure [T] is indicated by the light shading. Note that the 
abrupt kinks in the contours of constant rrij^^ are the result of level crossings where the neutralino 
goes from being over 95% wino-like (along y ~ 1) to over 95% bino-like (y ^ 1.4 and y ^ 0.6). The 
region with a wino-like LSP is indicated in Figure |2] by the nearly horizontal dashed lines. For the 
case of Mo = 500 GeV there is also some parameter space where the LSP is a mixture of Higgsino 
and wino, indicated by the hatched region. The analogous region for the case of Mo = 1000 GeV 



14 



would lie within the phenomenologically forbidden region. 

Some of the remaining superpartner masses are given in Figure [Sj again for the cases of Mq = 
500 GeV (left panel) and Mq = 1000 GeV (right panel). As before, the lightly shaded region is 
excluded. The upper panels show contours of constant gluino mass (solid heavy lines) and constant 
(lightest) chargino mass (dashed lines) for the entire {x, y} plane. The curved shaded region 
represents the area where 1.5 < mA/mj^_^ < 2.5. In this region we expected enhanced annihilation 
of relic neutralinos though the pseudoscalar Higgs resonance - particularly along the curved dashed 
contour for which rriA = 2m . Note that we have restricted our attention to the region in which 
both X and y lie between zero and one. This is the range explored in the mirage family of models. 
In the lower two panels we plot several superpartner masses as a function of the parameter x where 
we take Nm = 0, i.e. a true "mirage" model. The phenomenology in the lower panels roughly 
corresponds to following the line x + y = 1 in the upper two panels. Note that the correspondence 
between the plots, while very good, is not precise. Strictly speaking, the line x+y = 1 is determined 
by the requirement that R = 1. But from the definition in (19) it is apparent that there are two 
solutions for this constraint - one is the mirage limit where Nm = 0. If we fix the number of 
messengers at a value Nm (as in the case of the top panels in Figure |3]) then the line x + y = 1 
forces the second solution in which the quantity in braces in (19) vanishes. This is necessarily a 
solution for which ag does not vanish. Nevertheless, the phenomenology is similar to that of the 
mirage model near this line, so the reader can use the line x + y = 1 as a reasonable proxy for the 
mirage model limit. 

Finally, we wish to point out that much of the theoretically forbidden region is an artifact of the 



choice (30) for modular weights. It is clear from (26) and (27) that the choice of modular weights 
will have 0(1) effects on certain key parameters. This, in turn, affects the low-scale physical masses 
and the derived value of the ^-parameter once electroweak symmetry breaking is imposed. The 
amount of the {x, y} plane available is therefore heavily dependent on these model choices, as was 
pointed out in |24| [25| |26] . We therefore show the allowed parameter space in Figure |4] for the 
"model-independent" scenario in which we set fj, = tua = 1 TeV and set all scalar masses to the 
value of the gluino soft mass M3 or to 1 TeV, whichever is larger. As — > all gauginos get 
extremely massive for fixed Mq value. Thus, near the x = axis there is always a region where 
= 1 TeV, indicated by the dotted contour, where a scalar particle would be the LSP for our 
choice of scalar masses. As the gluino gets increasingly massive in the x ^ limit it induces large 
radiative corrections to the Higgs potential, producing regions where mh < 100 GeV (upper left 
corners) or where m\ < (lower left corners). The exact location of these contours depends on 
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Figure 3: Key Superpartner Masses for the "Model-Dependent Scenario." Upper panels give various 
masses for the {x, y} plane, while the lower panels plot masses along the line x + y = 1 with Nm = 0. Panels on 
the left take Mq — 500 GeV while those on the right take Mq = 1000 GeV. The lightly shaded region in all plots is 
phenomenologically forbidden. The darker shaded region in the upper plots shows the area where 1.5 < m^i/mjy^ < 
2.5, with tua = 2mjy^ given by the curved dashed line. In the upper plot heavy solid lines are contours of constant 
rrig while dashed lines are contours of constant m^^ . All masses are in GeV. 
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y Mo = 500 GeV y Mq = 1CX)0 GeV 




Figure 4: Allowed Parameter Space for the "Model-Independent Scenario." Gaugino mass contours are 
the same as Figure [l] The darker shaded region is the area in which the gluino is the LSP. The dotted line on the 
left of each plot is the contour where m~ = 1 TeV. The hatched region in the lower left of each plot indicates 
where m\ < and no EWSB occurs for a value fi = 1 TeV. We have indicated the area in the upper left where 
rrih < 100 GeV due to large radiative corrections. 



the choices made for fi and the soft scalar masses. 

As we have seen the allowed parameter space depends on the scalar sector of the theory, which 
is set by the modular weights assigned to the multiplets of the MSSM. The discrete choices for 
modular weights can be seen as selecting (a) the relative size of the gaugino masses versus scalar 
masses and (b) determining the degree of non-universality in the scalar masses between the Higgs 
sector and the matter sector. The former is akin to the choice of mo versus mi/2 in mSUGRA 
models, while the latter is related to certain non-universal extensions of these mSUGRA models. 
In order to capture this physics, in the next section we will consider the dark matter signals of the 
LSP neutralino in terms of the specific benchmark case of (30 ) over the space defined by {x, y, Mq}. 
But we will also occasionally relax this heavy constraint and randomize over all other parameters 
that determine the model. Together these two strategies will give a reasonable depiction of the 



dark matter implications of the gaugino masses in ( 22 ) . 
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3 Survey of Dark Matter Signatures 



3.1 Thermal Relic Density 

The phenomenology of relic neutralino dark matter consists of two largely disconnected aspects: 
the cosmological density of relic neutralinos and the physics of these neutralinos here in our galaxy 
today. These aspects are logically distinct, but can be related in any comprehensive theory of the 
origin and nature of dark matter. The physics of neutralino annihilation, or neutralino interaction 
with terrestrial dark matter detectors, depends on not much more than the properties of the LSP 
itself: its wave-function and mass. It also depends on its density in our local halo and certain other 
astrophysical properties independent of the nature of the LSP. The cosmological relic abundance, 
on the other hand, can be determined by this information only in certain special cases. More 
often it depends on the physical masses of other particles (supersymmetric states and Higgs states) 
that were part of the relativistic plasma at some earlier, hotter time in the cosmos. Crucially, 
the calculation of the cosmological relic density depends on certain assumptions about the thermal 
history of the universe. 

The most important inputs to this thermal relic abundance calculation are the mass of the LSP 
and the mass difference between it and the next lightest superpartner. The thermal relic abundance 



of the LSP neutralino for the theory with modular weights ( 30 ) is given in Figure [5] as a function of 
the dimensionless parameters x and y for the cases Mq = 500 GeV (left panel) and Mq = 1000 GeV 
(right panel). Here and throughout, all calculations are performed using the computer package 
DarkSUSY 5.0.4 [IS] after computing the physical mass spectrum throughout the parameter space 
using SuSpect 2.4. The WMAP three-year data [33 is best fit by a relic density in the range 

0.0855 < 17;^h2 < 0.1189, (31) 

at the 2a level. In our figures we will expand this region somewhat to a "WMAP preferred" region 

0.07 <n^h^ <0.U, (32) 



which will be slightly easier to resolve than the very narrow band in (31). 

In Figure [5] this "WMAP preferred" region is indicated by the narrow red shaded region. Inside 
the parameter space with a wino-like LSP, and within the ^-funnel region, the relic density drops 
below the value of ^x^"^ ~ 0-025, indicated by the yellow shading. The intermediate regime of 
0.025 < n^h.^ < 0.07 is indicated by blue shading. The green area m the plot has 0.14 < Q^h^ < 1 
while the gray area has ^x^'^ ^ ^- '^^^ transition from heavily wino-like LSP to heavily bino-like 
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Mo = 50O GeV y Mq = 1000 GeV 




Figure 5: Thermal Relic Density in Deflected Mirage Mediation. Left panel takes Mq = 500 GeV, right 
panel takes Mq = 1000 GeV. The "WMAP preferred" region of 0.07 < fi^h^ < 0.14 is here indicated by the narrow 
red shaded region. For the wino-like LSP and within the A-funnel region the relic density drops below the critical 
value of fi^h^ — 0.025, indicated by the yellow (very light) shading. The remaining regions are 0.025 < Q^h^ < 0.07 
(blue), 0.14 < n^h^ < 1 (green) and D,^h^ > 1 (gray). 



LSP occurs near y ~ 1.4 where the relic density rapidly changes from far too low to far too high. 
In between these extremes there exists a narrow region with n^h"^ ~ 0.10. 

The precision with which the cosmological abundance of cold dark matter can be inferred 
from the cosmic microwave background is remarkable. Yet despite the narrow window of (31) we 
cannot say that all of this cold dark matter is composed of neutralinos, nor that it was produced by 
standard thermal mechanisms. Thus we must be careful not to immediately exclude those parameter 
sets for which the thermal calculation yields ^^^^^ outside the range of (31). For example, an 
under-abundance of relic neutralinos may simply indicate the presence of non-thermal production 
mechanisms for stable LSPs |48| H21 ED] • Such mechanisms are especially well-motivated in string- 
inspired contexts such as in the KKLT model framework |5H |52] . An over-abundance of relic 
neutralinos is a greater cause for concern, but here too mechanisms exist which can bring the current 
relic density in line with the results from WMAP |53| [55] . We prefer to remain agnostic on the 
issue as this is not the physics of interest to us. In fact, the physics of relic neutralino detection can 
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Parameter 


Min Value 


Max Value 


Parameter 


Min Value 


Max Value 


Mo 


100 GeV 


2000 GeV 


A'mess 


10^ GeV 


10^^ GeV 


am 





5 




-5 


5 


tan/3 


2 


40 




1 


3 



Table 2: Parameter Scan Ranges for Deflected Mirage Mediation Models. In an effort to average over 
possible model-dependent effects we have produced 1000 deflected mirage mediation models by randomly choosing 
input parameters from within the ranges indicated. An additional random choice of modular weights for the matter 
and Higgs multiplets was also made to complete the model. 



be separated from the cosmological density of neutralinos to the extent that it is only the local halo 
density profile that is relevant. Of course the two will be related once a set of assumptions about 
the physics of the early universe are specified. We will, however, make one exception: for cases with 
n^h^ < 0.025 (indicated by the yeh ow shading in Figure 5^ it becomes difficult for the relic particle 
in question to account adequately for the local halo density of our galaxy [1]. Therefore when 
calculating observable quantities that depend on the relic neutralino number density n^^ present 
in our galaxy (or the energy density p-^ = my.n^) we will rescale the assumed local density of 
(p^)q = 0.3 GeV/cm^ by the multiplicative factor = Min(l, 



Displaying results in the {x, y} plane defined by (21) is useful for making contact with certain 
theoretical limits, but requires fixing a number of important input parameters. We will therefore 
also show results of a scan over parameter inputs for the defiected mirage mediation paradigm. 
The quantities W6 vary arc listed in Table |2| Note tliat tlie input parameters Aj-nggg and 771-3/2 are 



obtained from the definitions in (12) and (13), respectively. For the modular weights we randomly 
chose from the possibilities 

HQ = nu = riD = riL = riE = or 1/2 

riHu = = or 1/2 or 1 (33) 

which include the cases with the widest range of allowed points at the electroweak scale [26]. We 
generated 1000 points, all of which have proper electroweak symmetry breaking and satisfy the mass 



bounds imposed in Section 2.3 For the sake of comparison we also generated 1000 points with 
Um = Og = and vanishing modular weights. These points represent unified modulus-mediated 
models with a phenomenology similar to traditional dilaton-dominated models \56\ |57] . 

We plot the thermal relic density ^^^^^ ^ ^ function of the LSP mass for both the mSUGRA-like 
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Figure 6: Relic Density for Models from Table [2| Thermal relic density Q,^la^ is plotted as a function of 
the neutralino LSP mass. Dark circles are the ensemble of deflected mirage mediation models defined by Table [2] 
Larger lighter circles are an ensemble of unified modulus-mediated models with am — ag = 0. Red crosses are the 
DMM models which faU into the "WMAP Preferred" region 0.07 < f^xh^ < 0.14. Models in the shaded region below 
nh^i„ — 0.025 have the local halo density rescaled before any further observable is calculated. 



case Om = Og = and for the ensemble of deflected mirage mediation models defined by Table [2] 
in Figure [6| Dark circles are the ensemble of deflected mirage mediation models, while the larger, 
lighter circles are the cases with am = ag = 0. Red crosses are the DMM models which fall into 
the "WMAP preferred" region ([32]). Note that the WMAP-preferred area for the DMM model set 
tends to cluster around m^y^ ~ 0(1 TeV), as opposed to the unified models which require a much 
lower-mass LSP. We will designate this group of models with red crosses in the figures which follow. 
Finally, we note that those model in the shaded region below Jlh^j^ = 0.025 have the local halo 
density rescaled by the factor r-^ = Min(l, J7-^h^/0.025) before any further observable is calculated. 
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Ref. 


Experiment Name 


Fiducial Mass [kg] 


Exposure Time [yr] 


Rio [counts/(kg yr)] 


[63] 


XENONIO 


5.4 


0.16 


11.54 


[65] 


XENONIOO 


170 X 0.8 


1 


7.35 X 10-2 




LUX 


350 X 0.8 


3 


1.19 X 10-2 


m 


XENONIT 


1000 xO.8 


5 


2.50 X 10-3 


[62] 


CDMS II 


3.75 


0.29 


9.18 


m 


SuperCDMS (SNOlab) 


27 X 0.8 


3 


0.15 


m 


SuperCDMS (DUSEL) 


1140 X 0.8 


5 


2.19 X 10-3 



Table 3: Rate Estimates for Various Experiments. The minimum threshold rate 7?io necessary to produce 
10 events in a given experiment for the fiducial mass and exposure time given is tabulated in the final column. Note 
that for XenonlO and CDMS II we use the experimentally quoted fiducial masses. For all other (future) experiments 
we assume a fiducial mass equivalent to 80% of the nominal quoted target mass. 

3.2 Direct Detection 

Given the various theoretical considerations outlined above, it is clear that terrestrial direct detec- 
tion experiments would provide the firmest evidence for the existence of weakly-interacting massive 
dark matter. The interaction rate of relic neutralinos with target nuclei depends on the local den- 
sity of the neutralino in the local halo. This number is more tightly constrained by experimental 
observation [H EH] and fits to large scale structure simulations of galaxies similar to ours. The rate 
depends heavily on the properties of the neutralino itself, but it also requires the input of certain nu- 
clear matrix elements that must be inferred from experiment. The values of these matrix elements 
are uncertain - often to a surprisingly large degree. The resulting uncertainty in the interaction 
cross-section can be as much as 50% in some circumstances |59| [HO] . Though our interest here is in 
the study of broad correlations between the parameters of the generalized mediation scenario with 
the properties of the relic neutralino, we should keep these uncertainties in mind when discussing 
the results of any calculation. 

We will focus here exclusively on two types of detector design: cryogenic germanium bolome- 
ters and dual-phase liquid/gas xenon detectors. These types are currently operational in at least 
one experiment and producing results - and future enlargements are envisioned for each. We use 
DarkSUSY with default settings for all nuclear form factors to obtain the differential rate of inter- 
actions per unit recoil energy on germanium and xenon We do so over a range of recoil energies 

^We have corrected a normalization error in the recoil rate routines and ensured that all results are rescaled (when 
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relevant to the desired experiment. As in |61j, we perform the integration of these rates using two 
possible energy ranges 

Xe : 5 keV < ^recoil < 25 keV 

Ge : lOkeV < -Erecoii < lOOkeV. (34) 

Note that these rate calculations include both spin-independent and spin-dependent contributions. 

Bounds on these recoil rates have been set by current direct search experiments such as CDMS II 
[62] and XenonlO |63j, using exposures of 316.4 kg-days on xenon and 397.8 kg-days on germanium, 
respectively. Expected backgrounds for both experiments are low, typically on the order of (or less 
than) 10 expected background events per year of exposure. Background rates at future enlargements 
are expected to be even lower. We will therefore compute the necessary interaction rates over the 



energy ranges in (34) to produce 10 events for a given mass and exposure time in a number of 
current or planned experiments based on xenon or germanium targets. We summarize the detector 
details we will assume in Table [3j For the current experiments a limit can be placed by requiring 
R ^ ^10) the rate for which 10 events would have been produced at each experiment for the quoted 
exposure. These numbers are given in the final column of Table |3] and indicated in Figure [7] by the 
various shaded regions. We find that none of the model points give rates which should have been 
detected in the initial runs of XenonlO and CDMS II. 

For the future experiments listed in Table [3] we computed the value of Rio assuming a fiducial 
target mass equal to 80% of the quoted (nominal) mass and the exposure time listed. The results 
(in units of recoils per kg per year) for the {x, y} plane are plotted in Figure [7] for xenon and 
germanium targets, respectively. Results for our random model survey are given in Figure [8j The 
yellow (light shading) in all panels indicates parameter space that will be inaccessible to any of the 
future experiments in Table [3| At roughly the 100 kg- year level of exposure a fair amount of the 
parameter space will be detectable, particularly for low mass scales. This is indicated by the green 
(medium shading) in all panels, for which 11.54 > Rf^ > 0.0735 counts/kg-yr and 9.18 > Riq > 
0.15 counts/kg-yr for xenon and germanium, respectively. This is within reach of the XenonlOO 
experiment after about one year of operation, or the SuperCDMS SNOLab proposal after about 
three years of operation. In about one ton- year of exposure - equivalent to about three years of the 
LUX experiment - the red (dark) shaded region in the xenon panels with 0.0735 > Riq > 0.0119 
counts/kg-yr will be probed. Finally, after approximately five ton- years in xenon the blue shaded 
region with 0.0119 > Rf^ > 0.0025 counts/k g-yr can be explored, or the red shaded region with 



necessary) by the thermal relic density of LSPs. 



23 



Xenon Target 



Mo = lOOOGeV 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 X 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 X 

y Germanium Target Mq = 500 GeV y Germanium Target 

2.0 



1.5 




1.0 



0.5 



0.0 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 X 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 X 



Figure 7: Neutralino Recoil Rates on Xenon (top) and Germanium (bottom). Left panels set Mo = 
500 GeV, right panels set Mo = 1000 GeV. Phenomenologically allowed areas are enclosed by the heavy dashed 
lines. Colored shading indicates the reach of future direct detection experiments, as computed in Table [S] Yellow 
in all panels indicates parameter space that will be inaccessible to any of the future experiments in Table [3] For 
xenon targets, green indicates 11.54 > R^q > 0.0735, red indicates 0.0735 > R^q > 0.0119 and blue indicates 
0.0119 > -Rio^ > 0.0025 in recoils per kg-year. For germanium targets, green indicates 9.18 > R^q > 0.15 and red 
indicates 0.15 > i?5^o° > 0.00219 recoils per kg-year. 
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Figtire 8: Recoil Rates in Direct Detection Experiments for Models from Table [2] Reach of current 
and future experiments in terms of the 7?io values of Table [s] are given for xenon-based experiments (left panel) and 
germanium- based experiments (right panel). Dark circles are the ensemble of deflected mirage mediation models 
defined by Table [2] Larger lighter circles are an ensemble of unified modulus- mediated models with Om = Og = 0. 



Red crosses are the DMM models which fall into the "WMAP preferred" region (321. The xenon-based experiments 



are (a) XenonlO (reported), (b) XenonlOO (1 year), (c) LUX (3 year) and (d) XenonlT (5 year). The germanium- 
based experiments are (a) CDMS II (reported), (b) SuperCDMS (SNOLab - 3 year) and (c) SuperCDMS (DUSEL - 
5 year). 



0.15 > Rf^ > 0.00219 counts/k g-yr in germanium. Note that these regions are related to the 
wave-function of the neutralino LSP and to the presence of light Higgs states in the spectrum. In 
the ^4- funnel region where the relic density is well below the ^^^^^ ~ 0.025 threshold the interaction 
rates have been rescaled, thereby requiring a much larger exposure for detection. 

The results of the survey of randomized models are given in Figure [8} The current experimental 
limits are given by the horizontal lines labeled (a) in the figures. Though a few of our unified models 
would have been detected at these experiments, none of the deflected mirage models would have 
produced at least ten recoils. The five ton-year limit is indicated by the horizontal lines labeled (d) 
for xenon (with Riq = 2.50 x 10~^ counts/kg-yr) and (c) for germanium (with i?gf = 2.19 X 10"^ 
counts/kg-yr). The majority of the points which satisfy the WMAP preferred bound of (32) can 
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Exposure [km^ yr] 


^10 [counts/ (km^ yr)] 


0.2 


50 


0.5 


20 


1.5 


6.7 


10 


1.0 



Table 4: Muon Flux Estimates for IceCube. Flux of muons needed to produce 10 signal events at IceCube 
assuming various exposures. 

be probed at this level of exposure. 

3.3 Muons 

The presence of relic neutralinos can also be inferred from experiments which seek to detect the 
products of neutralino annihilation processes. The rate for neutralino pair annihilation will be quite 
low except in those areas where the present-day relic density is high. For indirect detection searches 
to be successful, therefore, they must be sensitive to objects that (a) are able to travel unimpeded 
for long distances from their point of origin, and (b) are not themselves copiously produced by 
mundane astrophysical processes. For the second point, it may be sufficient to demand that some 
property of the object (such as the differential flux as a function of energy) be significantly different 
from estimates of the astrophysical backgrounds. Any calculation that seeks to predict these signals 
will be faced with uncertain theoretical inputs - particularly the density of neutralinos at the extra- 
terrestrial location of interest. 

A good place to start, therefore, is at a location that is well-understood such as at the center of 
the sun or earth. Relic neutralinos can become gravitationally trapped in the core of the sun, sig- 
nificantly enhancing their probability for annihilation. Of the annihilation products, some fraction 
of the neutrinos can eventually exit the sun and be detected in experiments such as IceCube |67j 
via conversion of muon neutrinos into muons. To calculate this rate we integrate the differential 
flux of conversion muons from solar-born as well as earth-born neutrinos over the energy range 
50 GeV < Efj_ < 300 GeV, assuming an angular resolution of 3 degrees. The nominal target area 
for IceCube is 1 km^, but the effective area for detection of neutrinos via muon conversion is smaller 
and can be cast as a function of the muon energy |68l [69| [70] . To determine the prospects for ob- 
serving muons at IceCube we estimate various values of exposure in km^-years needed for 10 signal 
events, $io, as given in Table |4| We note that the estimated background of muon events in one 
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Figure 9: Upward Going Muon Rates at IceCube. Left panel shows the IceCube reach for Mq = 500 GeV. 

The red shaded region is detectable with 0.2 km'^-years of exposure. Green and blue shaded regions would require 0.5 
and 1.5 km^-years of exposure, respectively. The yellow region gives a visible signal after 10 km^-years of exposure. 
The gray region has a flux below 1 muon per km^-yr and is likely undetectable at IceCube. The right panel gives 
the muon flux at IceCube for the random models deflned by Table |2] Points above the horizontal line at — 400 
muons/km^-year would likely have given a signal in the 22-string IceCube data. The ultimate reach for IceCube is 
approximated by the lowest horizontal line at = 1 muon/km^-year. 

year of IceCube data taking is expected to be 0(10) |70] . 

In Figure [9] we show the number of observed muons at IceCube in units of counts per km^-yr, 
combining those from neutrahno annihilation in the sun with those in the center of the earth. The 
left panel shows the regions in the {x, y} plane satisfying the condition <I> > <I>io for the values 
listed in Table |4] for the case when Mq = 500 GeV. For higher values of Mq the (already low) muon 
rates drop rapidly and little of the parameter space would be visible at IceCube. This is evident 
in the right panel, which shows the muon rate in units of muons/km-^-year for the random model 
set of Table [2j Recent results using a 22-string detector at IceCube constrain the muon flux to be 
below 300-500 muons per km^-yr [7T] . None of the points in the left panel of Figure [o] produce 
a flux this large, though a few of the randomly-generated points would have produced a signal at 
IceCube. We conservatively place the limit at 400 muons per km^-yr, indicated by the horizontal 
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line in the right panel of Figure [9j 

For low values of the mass scale Mq there are some points which can produce a flux of 50 muons 
per km^-yr. These points are indicated in the left panel Figure [o] by the red shaded region. The 
three smaller flux values of Table [9] are indicated by the green, blue and yellow regions, respectively. 
This is roughly congruent with the parameter space that gives rise to a mixed Higgsino/gaugino 
LSP and a light supersymmetric Higgs sector. For the case where Mq = 1000 GeV and there is very 
little region with a mixed LSP we expect the flux of muons at IceCube to be well below 1 muon 
per km^-yr throughout most of the parameter space. The last two lines in Table |4] are indicated 
by the lower two horizontal lines in the right panel of Figure [9| Though the DMM models tend 
to give a larger flux of muons at IceCube than unified models with an equivalent LSP mass, the 
prospects for a visible signal at IceCube remain limited for most of the parameter space. 

3.4 Photons 

At distances further from the earth one can search for photons in the gamma ray energy regime, 
which travel largely unimpeded from their source. This allows gamma ray observatories to con- 
centrate on areas of the sky likely to have a high relic neutralino density - such as the center of 
our galaxy. Photons can be produced as part of decay chains or directly through loop-induced 
diagrams. The former contribute to a continuous general spectrum of photons and are thus more 
difficult to distinguish from astrophysical backgrounds. Direct production of photon pairs - or pro- 
duction of a single photon in association with a Z-boson - offers the possibility of a monochromatic 
spectrum that can more easily be distinguished from the background at the cost of a reduction in 
rate relative to the continuous photon rate. We will consider searches for both phenomena in this 
subsection. 

The calculation of the flux of gamma rays observed from the direction of the galactic center 
depends on the microscopic physics of the neutralino and its interactions, but also very strongly 
on the macroscopic physics of the halo profile assumed for the galaxy. The latter is conveniently 
summarized by a single parameter J (AO) where AO represents the solid angle resolution of the 
observatory. Here we will specifically consider the Fermi/ GL AST experiment for which we assume 
an angular resolution of AO = 10~^ sr. We will consider the commonly adopted profile of Navarro, 
Frank and White (NFW) \72^ which gives a value of J (AO) = 1.2644 x 10"^, as well as a modified 
NFW profile that takes into account the effect of baryons (a so-called "adiabatic compression" (AC) 
model) \73\ [71] for which we have J (AO) = 1.0237 x 10^. The large variation in calculated rates 
for different halo profile models should always be borne in mind when discussing "predictions" 
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Exposure [m^ yr] 


Halo Profile 


^iQO [counts/(cm^ s)] 


1 


NFW 


5.79 X 10"^° 


5 


NFW 


1.16 X 10~^° 


1 


NFW+AC 


7.15 X 10-12 


5 


NFW+AC 


1.43 X 10-^^ 



Table 5: Reach Estimates in the Continuous Gamma Ray Flux for the Fermi/GLAST Experiment. 

The quantity $ioo is the flux needed for 100 signal events at Fermi/GLAST assuming various exposure times and 
halo profiles. 

for gamma ray signals of any particular model. Nevertheless, we can study the dependence of 
the calculated rate on the underlying parameters and scale the normalization of that rate to any 
particular halo profile via the ratio of J {AO,) values, as needed. 

For the diffuse signals we use DarkSUSY to compute the differential photon flux from neutralino 
aimihilation in the halo in units of photons/cm^/s/GeV. We do so over the photon energy range 
relevant for the Fermi/GLAST experiment of 1 GeV < < 200 GeV, with a step size of 1 GeV. 
Note that the upper limit on the Fermi detector sensitivity is typically well below the mass of the 
neutralino LSP for much of the parameter space we are studying. We then numerically integrate 
these values to obtain an overall photon flux in units of photons/cm^/sec. Using the effective area 
of 1 m^ for the Fermi detector we compute the minimum flux $ioo necessary to produce an event 
rate of 100 signal photons from the galactic center in a given length of exposure. As an example, 
using the NFW proflle one obtains a value of $ioo = 5.79 x 10^^*^ photons/cm^/sec for one year of 
running, where we assume data collection for 200 days per year. If we require 100 signal photons 
over a five year mission this becomes $ioo = 1-16 x 10^^'^ photons/cm^/sec. To take into account 
the possibility of a more generous halo profile such as the NFW profile with adiabatic compression, 
we can simply rescale this reach by the ratio of J (AQ) for the two profiles. In this case this implies 
a reach to lower fiuxes by a factor of approximately 80. The corresponding $ioo values are listed 
in Table E 

Translating this into a discovery reach for Fermi/GLAST requires making an assumption about 
the background sources for gamma rays in the 1 GeV < < 200 GeV energy regime. We will 
use the estimate of [SH] which gives a background fiux of ^I'bkgrnd = 5.06 x 10^^*^ photons/cm^/sec, 
or 0(100) photons/year for the aperture and energy range of Fermi/GLAST. This coincides with 
sensitivity measures on the photon fiux at Fermi / GLAST often quoted in the literature \75\ [76] . In 
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Figure 10: Reach of the Fermi/GLAST Experiment for Integrated Gamma Ray Flux. Left panel is for 
Mo = 500 GeV, right panel is for Mo = 1000 GeV. Phenomenologically allowed areas are enclosed by the heavy 
dashed lines. For Mo = 500 GeV the green shaded region has 5.79 x 10"^° > 4>ioo > 1-16 x 10""^° photons/cm^/sec 
and therefore may be visible at Fermi/GLAST after five years of exposure with the NFW profile. If the fluxes are 
rescaled to the NFW+AC profile we obtain the reach shaded in blue (one year of data-taking) and yellow (five years 
of data-taking). The gray shaded area has an effective flux below — 0.014 x 10^^" and is unlikely to yield a visible 
signal for diffuse gamma rays at Fermi/GLAST. 



addition to looking for excesses in the overall rate, the Fermi experiment will also be able to use 
the shape of the observed spectrum to augment signal extraction. Furthermore, there is reason to 
believe that point sources near the galactic center can be identified and subtracted from the signal 
to further enhance signal significance [77] . We will therefore optimistically take the calculated flux 
<&ioo as a measure of the reach of the Fermi experiment in the deflected mirage mediation model 
space. 



The result of this reach analysis is given in Figure 10 for the cases of Mq = 500 GeV (left panel) 
and Mq = 1000 GeV (right panel). Detecting a gamma ray signal in the deflected mirage mediation 
model will be challenging, even with a relatively favorable set of assumptions about the halo profile. 
For the low mass case of Mq = 500 GeV there is a small region near the A-funnel (green shading) 
for which a signal is possible after five years of data-taking at Fermi/GLAST with the NFW profile. 
This coincides roughly with the region with a mixed Higgsino/gaugino LSP. No such region exists 
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Figure 1 1 : Integrated Gamma Ray Flux for Models in Table [2] Dark circles are the ensemble of deflected 
mirage mediation models defined by Table [2] Larger lighter circles are an ensemble of unified modulus-mediated 
models with am = Qg — 0. Red crosses are the DMM models which fall into the "WMAP preferred" region (32 I. 



Solid lines indicate Fermi/GLAST reach ($100) after one year and five years of exposure. The dashed lines indicate 
the same after rescaling all fluxes using the NFW + AC profile. 



in the higher mass Mq = 1000 GeV case. The blue and yellow regions in both panels are the reach 
at Fermi/GLAST after one year and five years, respectively, if the flux calculation was rescaled to 
the NFW + AC profile. Roughly speaking, with this set of halo assumptions it should be possible 
to detect a diffuse gamma ray signal over most of the parameter space involving a wino-like or 
mixed Higgsino / gaugino LSP. The gray area in the figures corresponds to the area with a bino-like 
LSP and typically has an effective flux below = 0.014 x 10"^'' photons/cm^/sec. These regions 
are unlikely to yield visible signals without additional favorable assumptions about the halo model. 
We note, however, that these rates would increase considerably if we were to assume a local halo 
density normalized to po = 0.3 GeV/cm'^ regardless of the calculated thermal relic density. 

The best prospects for diffuse gamma ray signals occur near regions that are ruled-out by 
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Exposure [m^ yr] 


Halo Profile 


<I>io [counts/(cm^ s)] 


1000 


NFW 


31.7 X 10"^^ 


100 


NFW+AC 


3.92 X 10"^^ 


500 


NFW+AC 


7.83 X 10-16 


1000 


NFW+AC 


3.92 X 10^16 



Table 6: Reach Estimates in the Monochromatic Gamma Ray Flux for a Generic ACT Experiment. 

Flux needed for 10 signal events at a generic ACT assuming various exposures and halo profiles. Note that we are 
here taking one year of data-taking to mean 365 days. 



certain scalars becoming the lightest supersymmetric particle. This, in turn, is a property of the 



modular weights assumed for the scalar sector of the theory. In Figure 11 we show the results for 
the integrated gamma ray flux in the range 1 GeV < E.y < 200 GeV for the random collection of 
models defined by Table [2] and modular weights ( 33 ) . Choosing different modular weights does not 
produce a large number of models with signals visible at Fermi/GLAST using the NFW profile. 
The one-year and five-year ^>ioo reach values are indicated by the solid horizontal lines. More 
encouragingly, the models that tended to satisfy the WMAP relic density constraints predicted 
values of ^>ioo that could potentially give observable signals at Fermi/GLAST if the NFW + AC 
profile were applicable. 

For the case of monochromatic signals, XX ~^ 77 oi' XX ~^ -^7; obtain the flux at a single 

2 

energy, = or = rriy. — respectively. As we have seen previously, for the deflected 
mirage mediation scenario the neutralino LSP can often be quite massive, particularly for the 
cases that satisfy the WMAP thermal relic density constraint. These monochromatic gamma ray 
signals are therefore likely to be at energies beyond the reach of the orbiting Fermi/GLAST satellite 
experiment. On the other hand, ground based atmospheric Cherenkov telescopes (ACTs) such as 
CAN CARGO [78], HESS [79], MAGIC [80] and VERITAS [81] have thresholds for photon detection 
in the 100 GeV range and can detect energetic photons up to 0(10 TeV). Though data-taking is 
restricted to dark, cloudless nights the effective area ^eff of the target detector is quite large. 

Unlike the case for continuum photons, the monochromatic signal has virtually no background, 
particularly at high photon energies. We therefore calculate the minimum flux <I>io to produce 10 
photons for a given exposure area x time, within the window defined by the energy resolution of 
a typical ACT. We take this energy resolution to be ge/E = 0.15. In Figure 12 we present the 
reach for a generic ACT with A^^ ~ 3 x 10^ cm^ using the estimates for <I>io given in Table km Here 
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Figure 12: Reach of the Fermi/GLAST Experiment for Monochromatic Gamma Ray Flux. Left panel 
is for Mo = 500 GeV, right panel is for Mq — 1000 GeV. Red shaded regions near the edge of the allowed area 
in the two plots have $10 > 31.7 x 10~^^ photons/cm^ /sec. Rescaled to the NFW+AC profile we obtain the reach 
shaded in green (100 m^-yr exposure), blue (500 m^-yr exposure) and yellow (1000 m^-yr exposure). Note that some 
of the indicated regions of the parameter space give photon energies below the lower threshold of — 100 GeV and 
are therefore invisible to our generic ACT. The gray shaded region has a monochromatic flux $10 < 0.39 x 10~^^ 
photons/cm^ /sec and is unlikely to give a visible signal without a more favorable set of halo assumptions. 



we have not included models with a neutralino LSP mass below 100 GeV, as the ACTs do not 
typically take data below this energy. As in Table |5] we convert the quantity <I>io into an estimate 
of the reach of our generic ACT over the model space for the NFW profile as well as for the NFW 
+ AC profile. Here we combine the photons from both the 77 and 7 Z signals as the 15% energy 
resolution will not be sufficient to resolve the two lines for LSP masses much in excess of 200 GeV. 

As in the case with the continuous gamma ray spectrum, detecting a signal from neutralino 
annihilation will be difficult without help from the halo profile. The red areas in the two plots 
have > 31.7 x 10~^^ photons/cm^ /sec and are potentially visible with the NFW profile given 
1000 m^-years of exposure. Note that the heavier mass cases (Mq = 1000 GeV) are easier to detect 
due to the sensitivity of ACTs to very energetic photons. If the monochromatic fiuxes are scaled 
to the values that one would obtain if we had used the NFW + AC profile then the area in green, 
blue and yellow shading (corresponding to 31.7 x 10^^^ > <I>io > 3.91 x 10~^^ photons/cm^ /sec. 
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Figure 13: Monochromatic Gamma Ray Flux for Models in Table [2] Dark circles are the ensemble of 
deflected mirage mediation models defined by Table [2] Larger lighter circles are an ensemble of unified modulus- 
mediated models with am = otg = 0. Red crosses are the DMM models which fall into the "WMAP preferred" 



region (321. The solid line indicates a generic ACT reach ($10) after 1000 m^-years of exposure. The dashed lines 
indicate the reach after rescaling all fluxes using the NFW + AC profile for 100, 500 and 1000 m^-years of exposure. 



3.91 X 10"^^ > $10 > 0.78 X 10"^^ photons/cm Vsec and 0.78 x 10"^^ > $10 > 0.39 x 10"^^ 
photons /cm^ /sec, respectively) would be visible after 1000 m^-years of exposure. This roughly 
corresponds to the area of wino and mixed Higgsino/gaugino LSP. 

In Figure 13 we show the results for the monochromatic gamma ray flux for the random col- 
lection of models defined by Table |2] and modular weights (33). The reach $10 for a generic ACT 
after 1000 m^-years of exposure is given by the solid horizontal line. Dashed horizontal lines give 
the reach after 100, 500 and 1000 m^-years of exposure, respectively, when all values of the flux 
are scaled up using the NFW + AC profile. The rates of annihilation through loop diagrams into 
mono-energetic photons are sensitive to the wave-function of the LSP. Hence the deflected mirage 



mediation models break into groups in Figure 13 determined by the makeup of the LSP. The unified 
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models, on the other hand, are nearly universal in having a bino-like LSP. We will return to this 
feature in our discussion in Section [4] As in Figure 11 we see that the WMAP preferred models 



are visible after a reasonable amount of ACT exposure provided the halo profile has a sufficiently 
high concentration at the galactic center. 

3.5 Anti-matter in Cosmic Rays 

Dark matter annihilations will also produce charged particles, of which positrons and anti-protons 
are the most promising for giving signals above astrophysical backgrounds. There has been exper- 
imental evidence for an excess of positrons in the 10-100 GeV range for some time |82| [83] . The 
recent data from the PAMELA experiment seems to confirm the previous results pi] . It is by now 
well appreciated that standard MSSM models can fit the signal from PAMELA only with some 
difficulty, and are even more challenged by the possible positron excess at even higher energies re- 
ported by the ATIC experiment [85j For that reason we will here only concern ourselves with the 
positron and anti-proton energy ranges relevant for PAMELA, roughly a few GeV up to 100 GeV. 

The excess over background was reported at higher energies, so we therefore consider the ten 
highest energy bins from the reported signal. We will take as the central value of these bins the 
values 

Ei = {5.55, 6.75, 8.25, 10.1, 13.05, 17.45, 23.85, 34.75, 53.2, 82.15} (35) 

for the energies (in GeV) of the electron or positron cosmic ray. Using DarkSUSY 5.0.4 (with all 
galactic propagation parameters set to their default values) we calculate the differential positron 
fiux from neutralino annihilation with the NFW halo profile in units of positrons/GeV/cm^/sec/sr 
and obtain an average fiux in each bin, denoted *^g+- We then form the ratio of positron fiux to 
positron and electron fiux as a function of the "boost factor" B via the relation 

MB) = —r-, ^ , (36) 

e+ e+ e~ 

where is the background positron fiux and (I)^"™^ are the contributions to the electron 

fiux from primary and secondary sources, respectively. Following |89| [90] we use analytic formulae 
to approximate these contributions from astrophysical processes 

$bknd ^ .37^ 

1 + 650^2.3 + 1500^4.2 y""' 



^As evidence for this statement one can consider the analysis performed in [861 1871 188] . among many other possible 



studies. 
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n 1 77—1.1 

= (38) 

1 + 11^0.9 + 3.2£;2.15 ^-^^^ 

= ^^if^ (39) 

where E represents the particle energy (normahzed to units of GeV). 

To quantify how weh a particular point in the parameter space of the DMM model can fit the 



PAMELA data, we use the function f{B) evaluated at each of the energy values in (35). The 
resulting quantities fi{B) are then compared to the observed PAMELA data points ff^^ using a 
chi-squared like function of the boost factor 

X^(i?)=E(^^^^^^^;^)'> (40) 

where the error values, Ui, are taken to be the average value of the full vertical error bar from the 
reported signal. We then determine the value of the boost factor B for which dx^/dB = and (40) 
is minimized. 

The results for the scan in the {x, y} plane for Mq = 500 GeV (left panels) and Mq = 1000 GeV 
(right panels) are shown in Figure [l4| The values of the boost factors that we obtain from mini- 
mizing x^{B) range from the relatively innocuous 0(10) to the preposterous 0(10^). The largest 
boost factors occur for the heavy bino-like LSP in the red shaded regions of the upper panels. The 
smallest boost factors are denoted by the yellow and purple shaded regions and occur for the wino- 
like LSP region (both panels) and the mixed Higgsino/neutralino LSP region (Mq = 500 GeV). 
Note that for these cases a large part of what is here termed a "boost factor" is in fact simply the 
result of a low thermal relic density. Invoking a non-thermal mechanism for explaining the present 
cosmological density of such wino-like LSPs would therefore result in 0(1) boost factors for many of 
the models in this part of the parameter space. This is consistent with previous results suggesting 
that a low mass wino-like LSP is the best candidate for an MSSM explanation for the PAMELA 
data pi] . 

The heavy bino-like LSP region fit the PAMELA data extremely well (purple and yellow shaded 
regions in the lower panels of Figure 14), but this is merely an artefact of the LSP being very 
massive - the PAMELA signal is then interpreted as a (massively-boosted) tail with a peak at a 
much higher energy level. More realistic is the wino-like LSP region. Though the boost factors were 
relatively small here, the calculated flux was not a particularly good fit to the shape of the resulting 
positron spectrum. Typical values of the minimized x^{B) were near unity (the blue shaded region 
in the lower panels of Figure 14), but the fi{B) values tended to fall with increasing energy like 
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Figure 14: Best-Fit Boost Factors and Minimized x^{B) Values for Fitting Anomalous Positron Data. 

Left panels set AIq = 500 GeV, right panels set Mo — 1000 GeV. Phenomenologically allowed areas are enclosed by 
the heavy dashed lines. In the upper panels the regions are separated according to best-fit boost values as follows: 
B > 10* red, 10** > B > 10*^ green, 10*^ > B > 10'' blue, lO" > B > 10^ yellow, and B < 10^ in purple. In 
the lower panels the regions are separated according to the minimized Xmin(^) values as follows: Xmin(-B) > 2 red, 
2 > xlUB) > 1.5 green, 1.5 > Xmin(-B) > 1 blue, 1 > xLniB) > 0.5 yellow, and xLniB) < 0.5 in purple. 
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the background, as opposed to showing an upward trend as seen by the PAMELA experiment. 
Nevertheless, when certain reasonable modifications to the cosmic ray propagation parameters in 
DarkSUSY are made, the shape of the spectrum can be made to fit the observations [02] • 

The PAMELA experiment has also reported data on anti-protons in cosmic rays which were 
consistent with expectations from background astrophysical sources as well as previous experimental 
data [93]. To compare the defiected mirage mediation parameter space to this data we compute 
the differential flux of anti-protons using DarkSUSY with the NEW profile and default diffusion 
parameters for the ten highest bins reported in [93], corresponding to anti-proton mean kinetic 
energies of 5.85, 6.98, 8.37, 10.1, 12.3, 15.3, 19.5, 25.9, 37.3 and 61.2 GeV. We use the estimations 
of Cirelli et al. [M] for the background flux from astrophysical processes. We consider a model point 
to be in conflict with the PAMELA data if the predicted signal flux in any one of the ten bins were 
twice as high as the reported experimental observation. None of the parameter space gives rise to 
an anti-proton flux of this magnitude. If, on the other hand, we were to naively apply the same 
boost factor to the anti-proton flux as we do to the positron flux then the majority of these models 
would be in conflict with the PAMELA anti-proton data. While anti-protons and positrons will 
propagate through the galaxy differently [95] , it is reasonable to assume that the fluxes will involve 
similar astrophysical backgrounds. Thus the two sets of data appear to be in conflict throughout 
the DMM parameter space unless modiflcations to the diffusion model from the DarkSUSY default 
values are made. 



3.6 Benchmarks 

Before concluding this section we will look back at the dark matter signatures for the speciflc 
benchmarks outlined in Table [T] in Section 2.3 To remind the reader, models A, B and C are 



specific examples presented in the original papers on deflected mirage mediation |29| |30] . Models 
D and E are chosen from the results of our parameter space survey, while model F is a point without 
gauge-charged messengers with am = 1 as in the KKLT model. The properties of the neutralino 
LSP, the masses of other key particles, and the values for the thermal relic density and other dark 
matter observables are given in Table [7j 

Properties of the LSP are given in the first block of entries in Table [7] The wave-function of 
the LSP is parameterized as 

Ni = NuB + NuW + NisH^ + N^H^, (41) 
which is normalized to N^i -\- + N^^ + = 1. The properties of the NLSP and key Higgs 
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Model A 


Model B 


Model C 


Model D 


Model E 


Model F 




1009 




626.6 




772.8 




257.2 




193.2 




710.8 




fB 


0.912 




0.989 




0.983 




0.914 




0.000 




0.932 




fw 


0.071 




0.001 




0.001 




0.030 




0.997 




0.004 




fH 


0.017 




0.009 




0.017 




0.057 




0.003 




0.064 




NLSP 


Ci 




Ci 








Ci 




Ci 




Tl 




m-NLSP 


1026 




743.8 




781.3 




302.6 




193.3 




762.8 




mh 


119.5 




120.6 




121.5 




112.1 




119.8 




119.9 




TUA 


1559 




1043 




1136 




473.5 




1518 




1136 






0.117 




0.891 




0.517 




0.088 




0.001 




0.537 




Rxe 


0.046 




0.022 




0.032 




1.182 




7 X 10"^ 




0.178 




RCe 


0.037 




0.018 




0.025 




0.881 




6 X 10"^ 




0.140 




*77 


1.778 X 10" 


17 


3.517 X 10" 


17 


2.734 X 10" 


17 


4.128 X 10" 


16 


2.071 X 10" 


15 


2.267 X 10" 


17 


%Z 


9.169 X 10" 


17 


3.244 X 10" 


18 


3.326 X 10" 


18 


8.868 X 10" 


16 


6.399 X 10" 


15 


1.157 X 10" 


17 


3>tot 


1.846 X 10" 


12 


2.260 X 10" 


12 


2.041 X 10" 


12 


1.449 X 10" 


10 


1.463 X 10" 


11 


4.319 X 10" 


12 




0.738 




0.194 




0.495 




11.766 




0.016 




3.616 






8.096 X 10" 


15 


2.155 X 10" 


14 


1.195 X 10~ 


14 


1.631 X 10" 


12 


2.824 X 10" 


14 


2.560 X 10" 


14 


*p 


7.438 X 10" 


11 


9.979 X 10" 


11 


8.923 X 10" 


11 


6.523 X 10" 


-9 


8.617 X 10" 


10 


1.887 X 10" 


10 


B 


4774 




5156 




5145 




102.9 




1147 




2510 





Table 7: Characteristics of Benchmark Models Prom Table [l] All masses are given in GeV and we denote 
wave-function components as fs = |A'"ii|^, fw ~ |A'^i2l^ and fn ~ \Nr.i\'^ + j7Vi4|^. Direct detection rates Ri have 
units of recoils/kg- year. Photon fluxes are given in units of photons/cm"^ /sec and "I>tot represents the diffuse gamma 
flux integrated from 1 to 200 GeV. The muon flux is given in units of muons/km^-year. is given in units of 
anti-deuterons/(GeV s cm"^ sr) and is computed at 0.25 GeV. is given in units of anti-protons/ (GeV s cm"^ sr) 
and is computed at 10 GeV. 



masses are given in the next block of entries. These can have a significant impact on the thermal 
relic density calculation, which is provided in the next group of entries. All of the benchmark 
models have a thermal relic density that is below unity, though only models A and D agree well 



with the WMAP results of (31). In computing the relic density values all the models receive some 
reduction in abundance due to coannihilation. This is particularly true for points A, C and E. For 
this last case the LSP is mostly wino-like and the computed relic density is well below the WMAP 
bound. It is also below the threshold value of Jl^h^jj^ = 0.025 and thus all subsequent dark matter 
signals have been rescaled to take this into account. Model D is an example of a point which resides 
in the A-funnel region where 2m~^ ~ rriA which also has effects on other dark matter observables. 
Apart from the wino-like model E, all of these cases would be probed by direct detection 
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experiments at the one-ton level, and most would give a sizable signal for detectors in the 100-300 
kilogram range. As previously mentioned, one should keep in mind that these rate values do come 
with large theoretical uncertainties associated with the nuclear matrix element factors which go 
into the calculation. Models D and F have the largest Higgsino content and hence the largest direct 
detection rates, with model D giving the most sizable rate due to its low mass LSP. We point out 
that the very low rate for model E is in part due to the rescaling of the local halo density we 
performed. If we had simply assumed a local halo density normalized to 0.3 GeV/cm^ these rates 
would be increased by a factor of 25. 

Gamma ray signals for the benchmarks are shown using the NFW halo profile. For the con- 
tinuous gamma ray signal we estimate that a flux on the order of 1-5x10^^*^ photons/cm^ /sec in 
the Fermi/GLAST photon energy range of 1 to 200 GeV is needed to observe a signal over the 
astrophysical background. Only model D gives rise to such a large flux of gamma rays. Each of the 
others, however, would be observable above background if the NFW profile with adiabatic compres- 
sion were employed. The wino-like LSP of model E gives a large signal for the two monochromatic 
gamma ray signals. Such a (combined) signal is just above the threshold for most ACT experi- 
ments and would give a visible signal with the NFW profile for 5000 m^-years of exposure - not an 
unreasonable amount given the size of most atmospheric Cherenkov detectors. For the NFW -|- AC 
profile it is likely that only models A, D and E will give rise to a detectable monochromatic gamma 
ray signal. Expectations for observation of these benchmarks at IceCube via the flux of muons 
from neutralino annihilation in the sun are also a bit mixed. Models D and F have the highest 
muon fluxes and ought to produce muons at rates above background for reasonable estimates of 
the effective area of the IceCube detector at these energies. The remaining models are unlikely to 
produce a signal even after 10 km-^-years of data-taking at IceCube. 

Anti- matter prospects fare slightly worse across our benchmark models. In all cases the naive 
estimate of anti-proton fluxes is well below the background expectation of 0(10^^) p/(GeV s cm'^ 
sr) at our reference point of 10 GeV of kinetic energy for the anti-proton. If one wishes to explain 
the PAMELA positron data using these models, however, one is confronted with rather large boost 
factors for all models save model D. Naively applying such boost factors to the anti-proton flux 
would place all of these models in tension with the anti-proton data, giving a signal flux of roughly 
an order of magnitude above the background. For completeness we have also included the flux of 
anti-deuterons computed at a reference point of 0.25 GeV for the D kinetic energy which is relevant 
for the planned GAPS experiment [96] . Only model D again is likely to give a sizable signal in this 
indirect detection channel. 
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4 Mirage Mediation versus Deflection Mirage Mediation 



In the previous section we looked at how a number of astroparticle/cosmological signals depend on 
the parameters of the deflected mirage mediation paradigm. We did this in terms of the theory 
motivated {x, y} parameter space where we saw that the detection prospects are good for general 
values of the parameter x when 0.4 ^ y ^ 1.4 for the modular weight choice (30 1 . We also averaged 
over the other more model-dependent parameters in the theory, presenting our results in the form of 
scatter plots, using a unified mSUGRA-like theory as a foil with which to compare the predictions 
of the DMM model. That DMM and unified models give different "footprints" should not be a 
surprise, since the importance of non-universality in gaugino masses for dark matter observables has 
been emphasized for some time |97l [98] . More interesting is the question of whether observations 
in the dark matter arena are capable of distinguished DMM from the (un-deflected) mirage pattern 
which preceded it. 

Looking back at the figures from the previous section we might conclude that the answer is 
certainly affirmative. After all, when we move away from the mirage prediction x + y = 1 we see 
strong variation in both the mass and composition of the LSP, resulting in dramatic changes in 
the thermal relic density of the neutralino. Moving away from the x + y = 1 limit can increase 
recoil rates at direct detection experiments and change the strength of signals at indirect detection 
experiments. But this is somewhat deceptive, since in those figures a number of important param- 
eters - particularly the overall scale Mq ~ were held constant. When we average over these values 



the distinction is far less clear. To exemplify this, consider Figure 15 in which we look at two of 
the more commonly considered observables in the literature: recoil rates on xenon targets and the 
integrated photon flux above 1 GeV in energy from the direction of the galactic center. In this 
figure we have fixed the value of Mq = 500 GeV. The dark points continue to represent a random 
collection of 1000 deflected mirage mediation models with parameters taken from the ranges given 
in Table [2} The lighter shaded points are now no longer a unifled model, but instead represent an 
additional 1000 points with the same ranges of Table [2] except that we flx Nm = 0. These are, 
therefore, mirage models with a similar overall mass scale. Though the range in LSP mass is more 
restricted in the un-deflected case, the signature range expected is nearly identical. These exam- 
ples could be replaced with many of the signals we considered in Section [3| Such "inclusive" style 
measurements in the cosmological arena will not in themselves be able to determine the presence 
of the gauge mediation component in the underlying theory. 

The reason for this degeneracy is readily found by considering another way to parameterize the 
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Direct Detection: Xenon Target 
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Figure 15: Comparison of Deflected Mirage Mediation with the Case of No Messengers. A random 
sample of models based on the ranges of Table [2] with Mo fixed at 500 GeV for the deflected mirage case (dark points) 
versus the same ranges but with Nm = (light points). The left plot gives the recoil rate per kg-year on xenon 
targets with sensitivities for XenonlOO after one year (solid line), LUX after three years (dashed line), and XenonlT 
after flve years (dotted line). The right figure gives the photon flux for all energies above 1 GeV from the galactic 
center. The solid line is the reach of Fermi/GLAST after flve years with the NFW proflle while the dashed line is 
the flve year reach with the NFW + AC proflle. 



gaugino masses in the DMM model. We can translate each of the instances of the running gauge 
coupling 5a (/u) in (11) to the high scale /Uuv = /"gut via the relation 
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The gaugino masses can now be written 
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(43) 



where the first term in the square brackets is precisely the case of gaugino mass unification at the 
high scale, as in the minimal supergravity paradigm. As in |2H l35] we see that departures from 
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Parameters 





a 


Mo [GeV] 


Mo [GeV] 




Nm X Ug 


Mirage 


6.52 


587 


1976 


1.5 





DMM 


6.52 


599 


889 


3.4 


-4.3 


Physical Quantities 




H% 


[GeV] 




<[cm2] 




Mirage 


98.1% 


1247 


0.13 


0.63 X lO-'^s 


6.51 X 10^12 


DMM 


99.5% 


1040 


0.10 


2.86 X 10-^5 


2.41 X 10-12 



Table 8: Two Degenerate Models from the Mirage Family. The mirage model and the DMM model have 
nearly identical values for the gaugino mass-determining parameters of equation (44 1. Both have nearly Higgsino-like 
LSPs of approximately a TeV in mass. Both give rise to a thermal relic abundance of LSP neutralinos in line with 



the WMAP preferred region of (32 I and have similar predictions for the observables in Figure 15 



universality are governed by the parameter am, but the slopes of these departure trajectories are 
altered relative to the case of simple mirage mediation. We can make this more explicit by defining 
new variables 
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such that ( 43 ) becomes 
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From (45) it is immediately apparent that despite the added complexity of the gauge-charged 
messenger sector the gaugino masses continue to depend on only a single dimensionless parameter 
in addition to an overall mass scale, here given by Moj^ In other words, for every mirage model with 
Nm = there exists a family of deflected mirage models with the same set of soft supersymmetry 
breaking gaugino masses at the electroweak scale. Any differences in phenomenology will therefore 
be due to the differences in scalar masses and the resultant EWSB parameters. 

The parametrization in (45) makes it easy to construct such look-alike "degenerate" points. 
Consider the two example points in Table [s] The (un-deflected) mirage model and the DMM 
model have nearly identical values for the key parameters a and Mo of (44) despite having very 
different underlying physics inputs. The mirage model has modular weights n = 1/2 for all fields 



Note that this variable is equivalent to Ms defined in ,22^ and utilized by [231 126 
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and tan (3 = 2.5 while the deflected mirage model has vanishing modular weights and tan /? = 27.8. 



Both are in accordance with the WMAP "preferred" relic abundance in equation ( 32 ) and have 



similar predictions for the two observables plotted in Figure 15, In particular, both would give an 



observable signal in one ton-year of exposure on a xenon target but neither would be detectable 
above background at Fermi/ GL AST with the NFW halo profile (though both would be very nearly 
detectable if the NFW + AC profile was used). Their predictions for other observables are also 
similar in magnitude. 

Though many such degenerate pairs can be constructed, it is not the case that there are no 
signals associated with relic neutralinos that are sensitive to the presence of gauge mediation in 
the underlying model. We will here merely point out two such examples. To do so we will change 
variables once again to a dimensionless quantity that is directly proportional to the number of 



gauge-charged messengers in the theory. Writing out (17) in terms of the explicit mass scales in 
the theory as 
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we may apply the definitions (12) and (13) to write (46) as 
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where we have introduced a new dimensionless parameter a' defined as 
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(48) 



Although specifying the value of o! does not uniquely determine the underlying model, all mirage 
models will have q' = 0, which makes this a convenient parametrization for our purpose here. When 
one plots inclusive counting signatures as a function of a' - such as the number of recoils per kg- 
year on xenon targets or the total integrated fiux of photons coming from the galactic center - the 



results are much as in Figure 15, But certain exclusive measurements reveal some correlation. For 
example, one can take advantage of the Fermi/GLAST experiment to measure the spectral profile 
of detected gamma rays. The shape of this spectrum carries information about the mass and wave- 
function content of the annihilating neutrahno LSP. We can attempt to take this information into 
account by using DarkSUSY to compute the differential photon fiux in 1 GeV increments over the 
energy range 1 GeV < < 200 GeV. From this information we create an interpolating function 
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Figure 16: Hardness of the Gamma Ray Spectrum from the Galactic Center as a Function of a'. Flux 
of photons in the energy range 60 GeV < -E7 < 200 GeV divided by the flux in the range 1 GeV < -B^ < 60 GeV is 
plotted for the signal photons versus the parameter a' of equation (481. The same ratio for the background is given 
by the heavy horizontal line. Note the logarithmic scale of the vertical axis. 



which is then integrated over six energy bins: 1-10 GeV, 10 - 30 GeV, 30 - 60 GeV, 60 - 100 GeV, 
100 - 150 GeV, and 150 - 200 GeV. 



In Figure 16 we take the ratio of the total flux over the three highest energy bins to the total 
flux over the three lowest energy bins and plot it as a function of the parameter a'. Over the range 
— 10 < q' < 10 (which is the regime of moderate ag values) the hardness of the photon spectrum is 
reasonably well correlated with the value of a' ~ x Og. Note that over this range of parameters 
we generally expect the gamma rays from neutralino annihilation to have a harder spectrum than 
those coming from background astrophysics sources, as indicated by the heavy horizontal line. The 



collection of models in Figure 16 was generated using a fixed mass scale of Mq = 500 GeV. The 
degree of correlation increases if we choose higher values for this parameter. 

More generally, the presence of the gauge-charged messengers allows the theory space to explore 
neutralino LSPs with differing wave-function composition from that of the mirage case with = 0. 



45 



100 



XI 

c 
o 
u 



10 



01 
Q. 

n 

E 

u 

e 
5 

a 0.1 

X 

^''o.oi 



Wlno-llke LSP * 



ACT(30D 571^ -yr) 



ACT <500 n^-yrj 



Hlggslno- like LSP 



...„.„-.. .:-jL.,-„jJ!.Jfa 

.i ' " 1 



BIno- like LSP 



100 



150 200 

Mji, (GeV) 



300 



500 



700 



Figure 17: Combined Monochromatic Line Flux as a Function of Neutralino Mass. The solid line 
represents the reach at a generic ACT for 500 m^-years of exposure with an NFW profile. The dashed line is the 
reach for the same exposure assuming an NFW + AC profile. 



This is already apparent from Figure 13 Such effects are known to have a large impact in the rate 



for loop-induced annihilation into monochromatic photons. In Figure 17 we give the combined flux 
at both the 77 and 7 Z energies for our random sample of DMM models with Mq = 500 GeV 
(dark points) and our random sample of mirage model (light points) at the same Mq value. When 
the overall mass scale Mq is held fixed the points tend to cluster as a function of the neutralino 
wave-function composition. In the mirage mediation model without gauge messengers the majority 
of models have a bino-like LSP with a smaller grouping of Higgsino-like LSPs. In fact, the maximal 
value of the wino-component of the LSP was typically 20%, and these occurred mostly in cases 
with vanishing modular weights for the matter multiplets and n = 1 for the Higgs multiplets. By 
contrast the DMM models are able to achieve a wino-like LSP across all modular weight choices. 
Assuming an NFW halo profile only those cases where the LSP is nearly 100% wino-like have a 
chance to be observed in a typical ACT experiment. Assuming an adiabatically-compressed profile 
allows most of the wino-like and Higgsino-like cases to be detected in 500 m^-years of exposure. 
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5 Conclusion 



Deflected mirage mediation, like the simpler mirage pattern which preceded it, are motivated 
by realistic string compactifications which stabilize all geometrical moduli while simultaneously 
addressing the vacuum energy problem of supersymmetry breaking. As such they are important 
cases for phenomenological investigation in their own right and have justly received much attention. 
More generally, if we are indeed on the eve of a new supersymmetric data era then it is of the utmost 
importance to begin the study of how this wealth of data can be used to understand the underlying 
principles of supersymmetry breaking. The deflected mirage paradigm provides a framework that 
is remarkably rich, allowing a smooth interpolation between supergravity (modulus) mediation, 
anomaly mediation, and gauge mediation. It is an important empirical question to ask how well 
we can isolate these various contributions to the transmission of supersymmetry breaking to the 
observable sector if all three are present at roughly the same mass scale. 

Ultimately we would like to answer this question in as model-independent a manner as possible, 
a point emphasized by Choi and Nilles [21] • As those authors point out, gaugino masses are a much 
cleaner window on the nature of supersymmetry breaking and transmission than the rest of the soft 
supersymmetry-breaking Lagrangian. This remains true in the deflected mirage paradigm as well. 
For this reason we have chosen to begin answering the question of the previous paragraph with 
dark matter observations which depend strongly on the gaugino sector and relatively less strongly 
on the scalar sector. Our results indicate that loosening the restrictions imposed by the original 
mirage pattern allows for an expanded parameter space in which the LSP can be predominantly 
wino-like, or mixed Higgsino/gaugino. This regime roughly corresponds to circumstances in which 
^ Nm 'X Oig ^ 10, or where A^mAmess = cm3/2 with c = 0(1). In such cases a number of direct and 
indirect detection experiments running now or scheduled to begin data taking in the near future 
have excellent prospects for detecting the relic neutralino. This is also the area in which a wino-like 
LSP with roughly the right properties to explain the PAMELA positron data can be found. For 
this statement to be borne out, however, it is still necessary to make modifications to the standard 
diffusion parameters and to imagine non-thermal production mechanisms to explain the present 
number density. It would be worthwhile to further investigate the utility of DMM-like models at 
fitting the PAMELA positron data. If the data does indeed seem to point toward a wino-like LSP, 
it would be interesting to see how a mixed-mediation wino arising from the DMM scenario would 
fare when a more sophisticated analysis is performed. 

Though the predictions from the deflected mirage mediation model depend strongly on the pres- 
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ence and magnitude of the gauge-mediation sector, it is not necessarily the case that observations 
in the dark matter arena are in and of themselves sufficiently powerful to determine this fact. That 
dark matter observations should suffer from such an "inverse problem" is perhaps not surprising. 
In fact, the gaugino sector of the DMM scenario is still a two-parameter family of models and is 
thus not fully general. Any trio of gaugino masses in the DMM theory space can be mapped to an 
equivalent trio on an un-deflected model. Thus dark matter signals will distinguish between the 
two paradigms only via the scalar sector, which serves to determine the LSP properties only via 
the clcctrowcak symmetry-breaking conditions. Using this fact to detect the presence of the gauge 
messenger sector will prove difficult, but some handle can be obtained via gamma ray signals. To 
make optimal use of this data we would prefer to have some knowledge of the relevant mass scales 
for the gaugino sector. For this it is crucial that the collider signatures for the deflected mirage 
mediation model be computed as a function of the contribution of the three mediation mechanisms. 
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